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THERMAL CORRECTION OF DEFORMATIONS IN A TELESCOPE MIRROR* 


By Marvin D. Rhodes 
Langley Research Center 

SUMMARY 

Orbiting astronomical observatories have the potential for making observations far 
superior to those from earth-based mirrors. In order for this performance to be real- 
ized, the contour of the primary mirror must be very accurately controlled. A prelim- 
inary investigation of the use of thermally induced elastic strains for correcting axi sym- 
metric deformations in space telescope mirrors has been presented. The relation between 
axial deformation and thermal inputs was determined by a finite -difference solution of the 
equations for thin elastic shells. 

The use of this technique was demonstrated analytically on a beryllium paraboloid. 
This mirror had 10 equally spaced thermal inputs and results are presented which show 
the nature of the temperature distribution required to correct deformations due to an 
acceleration-type loading. 

INTRODUCTION 

Orbiting astronomical observatories have two potential advantages over earth -based 
telescopes (ref. 1): first, the entire electromagnetic spectrum is available for observa- 
tion of stellar sources, and second, orbiting astronomical observatories will not be affected 
by air turbulence which limits the resolution of earth-based telescopes. Away from this 
turbulence the resolution of astronomical observations will be increased by at least 1 order 
of magnitude, which will permit telescope operation near the diffraction limit (ref. 2). 

Several studies of space telescopes have been conducted and all have recommended 
the use of a Cassegrainian optical system (refs. 3 and 4). These studies indicated that the 
primary mirror is the critical element in this system and deviations from the design shape 
must be controlled to extremely close tolerances. If close tolerances cannot be main- 
tained, the resolution of space telescopes may be no better than earth-based observatories 
and their usefulness will be severely restricted. 


*The information presented herein was included in a thesis entitled "Thermal Active 
Optics Technique for Correcting Symmetrical Distortions in Space Telescope Mirrors" 
submitted in partial fulfillment of the requirements for the degree of Master of Science in 
Mechanical Engineering, University of Virginia, Charlottesville, Virginia, September 1970. 


Two types of primary mirrors have been proposed for use in space telescopes 
(ref. 5). The first type is a passive mirror which would be designed to retain the proper 
contour without correction for the life of the telescope. This design may not be feasible 
since intolerable deformations may be caused by the relaxation of residual stresses. 

Also, actual launch and space environmental loads on the mirror are difficult to predict. 
The second type of primary mirror is known as an active optics system and uses a rela- 
tively thin mirror which is permitted to deform moderately under operational loads. A 
figure error sensor would monitor the deformations and activate a closed -loop control 
system to apply corrective loads to remove the deformations. 

One active optics system which uses an array of force actuators in conjunction with 
a closed-loop control system to provide the corrective loads has been investigated analyt- 
ically. This system has been experimentally verified (ref. 6) by using a thin, deformable, 
spherical mirror 76.2 cm in diameter. Another type of active control system that has 
been considered uses elastic deformations introduced by differential heating to maintain 
the proper contour. A limited thermal system has been successfully used on the earth- 
based telescope mirror of the Observatoire de Haute Provence (ref. 7). In the present 
paper a computer numerical analysis for predicting axisymmetric, thermally induced, 
corrective deformations in a primary mirror surface is presented. This analysis makes 
use of stored influence coefficients of deformations due to unit thermal inputs at various 
locations. The influence coefficients used herein were obtained from a finite -difference, 
linear shell analysis. The present analysis is demonstrated on a sample mirror configu- 
ration and the merits of two patterns of thermal inputs are investigated. In addition, the 
temperature distribution necessary to correct deformations from a constant acceleration 
loading is presented. 

SYMBOLS 

The units used for physical quantities defined in this paper are given in the Inter- 
national System of Units (SI). Factors relating this system with the U.S. Customary Units 
are presented in reference 8. 

A amplitude of applied temperature distribution 

c coefficients of flexibility matrix 

D diameter of telescope mirror 

f focal length of paraboloidal telescope mirror 
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k 


coefficients of stiffness matrix 


R principal radius of curvature of paraboloidal shell middle surface 

r radial coordinate 

r m radius of telescope mirror 

T applied temperature distribution at reference surface of mirror 

a coefficient of linear thermal expansion 

A meridional distance, see figure 4 

6 deformation of middle surface of shell in axial direction 

6a allowable deformation, 0.05 /im 

£ coordinate measured along shell meridian 

co principal curvature of middle surface of paraboloidal shell 

Subscripts: 

ij deformation at point i due to a unit load at point j 

n maximum number of control points 

6 circumferential direction 

£ meridional direction 

GENERAL CHARACTERISTICS OF THERMAL ACTIVE OPTICS SYSTEMS 

A sketch of a typical earth orbiting space telescope is shown in figure 1. This 
model is similar to the type discussed in reference 3. The basic configuration consists 
of two cylindrical shells which are attached to the telescope cabin. The cylindrical 
shells enclose the main optical elements which are the primary and secondary mirrors. 
Attached to the outer shell is a system of doors that prevents sunlight from falling on the 
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Station 

docking 

mechanism 



Figure 1.- Sketch of an earth orbiting space telescope. 

optical system during maneuvers. The inner shell is a thermal shield to reduce solar 
heating loads on the primary mirror. All optical imaging devices and sensing instru- 
ments are contained in the telescope cabin. The cabin provides the necessary environ- 
ment for manned support. 

A schematic of a thermal active optics closed-loop control system is shown in fig- 
ure 2. A complete description of a similar system utilizing force actuators can be found 



Figure 2.- Schematic of thermal active optics closed-loop control system. 
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in reference 6. In the present report, the primary telescope mirror is a thin parabolic 
shell of revolution. Deformations in this mirror are detected by a figure error sensor 
by using interferometric techniques. An analyzer then calculates the instantaneous 
amplitude of the temperature at the control stations necessary to correct the deforma- 
tions. An acceptable level of distortion will depend upon the wavelength of the light being 
viewed and upon the desired resolution. 

In general, thermal loads, that is, controlled temperature inputs, may be applied by 
conduction pads attached to the rear surface or by circulating fluid through channels in 
the interior of the mirror. For proper control and distribution, both heating and cooling 
loads may be necessary. Therefore, such a control system would probably be used in 
conjunction with a heat exchanger system. To correct asymmetric deformations the area 
over which the thermal load is applied will vary with the radial and circumferential loca- 
tion of the control station. The area and location of actuators for discrete control of lin- 
ear systems has been previously investigated and is reported in reference 9. The loca- 
tion of error sensors has also been examined in reference 9. 

SAMPLE MIRROR CONFIGURATION AND LOADING 
Mirror Geometry and Properties 

In order to demonstrate the use of a thermal optics control system an analytical 
model was established. A sketch of the mirror used in this analysis is shown in figure 3. 
The mirror is a thin paraboloidal shell having a focal -length —diameter ratio (f/D) of 4:1. 
This value is within the range of f/D currently being considered for space telescopes. 
Since no firm design of the telescope has been formulated, certain assumptions concerning 



Figure 3*- Sketch of cross section of paraboloidal 
telescope mirror. 
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the mirror geometry were necessary. For example, the diameter of the central hole will 
depend upon the telescope optical system and manufacturing consideration. The hole was 
assumed to have a diameter equal to 1/10 the diameter of the mirror. In the weightless 
environment of space, only a very thin reflective surface would be necessary for the pri- 
mary mirror. However, practical considerations of manufacture and launch will require 
the mirror to be sufficiently thick to withstand grinding, polishing, and handling in a grav- 
ity environment. For this investigation the thickness of the mirror was 1/100 of the 
diameter (ref. 5). The thickness was assumed to be constant in both circumferential and 
meridional directions. 

Metal mirrors are ideally suited for thin one-piece construction because they have 
high stiffness -weight ratios. Beryllium has one of the highest stiffness-weight ratios of 
any structural metal and is currently being considered as a prime candidate for telescope 
mirror construction (ref. 10). The material properties and dimensions for this mirror 
are shown in table I. This material was considered to be both homogeneous and isotropic 
in this analysis. 


TABLE I.- MATERIAL PROPERTIES AND DIMENSIONS 
OF BERYLLIUM MIRROR 


Property Value 

Modulus of elasticity, GN/m^ 276 

Poisson's ratio . 0.08 

Coefficient of thermal expansion, a , K"1 ............ 1.24 X 10“^ 

Diameter, D, m ................... . 3.0 

Focal length, f, m . 12.0 

Thickness, cm ................ 3.0 


The mirror was restrained in the axial direction and unrestrained in the radial 
direction at the outer edge (fig. 3). This support system is similar to the three -point 
tangent bar mounting suspension considered in reference 3. The central cutout portion 
of the mirror was unrestrained. 

Assumed Thermal-Active Control System 

For illustrative purposes, only steady-state symmetrical deformations of the pri- 
mary mirror were considered. As indicated previously, an acceptable level of deflection 
will depend upon the wavelength of the light source being viewed and the desired resolu- 
tion. For this analysis a maximum allowable deformation of 0. 1 wavelength of light was 
selected (ref. 3). 
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Ten points along the meridian of the mirror were selected as the control stations 
for the thermal loads. These thermal loads are rotationally symmetric about the shell 
axis. The use of more control stations would only have generated a larger matrix of 
influence coefficients and would have added little to the demonstration of the technique. 

For an actual control system a larger number of stations may be desirable to increase 
the mirror quality. Sources with a controlled temperature distribution were selected to 
provide the thermal inputs. The controlled sources were applied to the reference surface 
of the mirror. In order to simplify the problem, only steady -state deformations were 
examined. 

Two control patterns were considered. In one pattern the control stations were 
equally spaced along the meridian and in the second they were located in equal annular 
areas. 

The controlled temperature sources applied to the reference surface of the mirror 
were assumed to have the form 



where i = 1, 2, . . ., 10. A sketch of these sources applied at a typical control station 
on the mirror surface is shown in figure 4. Sources of this form were chosen because 
they represent a compromise between a realistic thermal load and one that is simple 
enough to demonstrate the analytical approach. The application of this thermal input will 
result in two-dimensional heat flow within the mirror interior. However, a preliminary 
examination has indicated that the interior temperature will not be significantly different 
from the applied back surface distribution. 



Figure 4.- Temperature distribution applied at control points along 
back surface of mirror for equal meridional spacing. 


7 



METHOD OF ANALYSIS 


In order to determine the influence of the thermal inputs it was assumed that defor- 
mations in the primary mirror are small and that linear shell analysis is applicable. The 
use of linear shell theory permits the application of the principle of superposition. Con- 
sequently the deformation at the stations may be expressed as a function of the applied 
thermal loads by the following linear equations: 


5 1 = C 11 T 1 +c 12 t 2 +■• • • c ln T n 
6 2 = c 21 T 1 + c 22 t 2 + • • • c 2n T n 


( 2 ) 


§n - c nlTl + c n 2 T 2 + . , . c nn T n J 

The influence coefficients cy specify the contribution of each temperature distribu- 
tion Tj to the deflection 6j at the control point 1. These equations may be written 
in the more convenient matrix notation as 


H =[c]|t| (3) 

The square matrix [c] is the flexibility matrix and the component terms are deflection 
influence coefficients. 

The thermal loads necessary to correct mirror deflections can be expressed by 
inverting the flexibility matrix and multiplying by the measured deformations. Thus 

|t| = [k]|s| (4) 

where 

frMc ]- 1 

and [k] is a thermal stiffness matrix. In order for the active optics system to perform 
properly it is necessary to determine accurately the coefficient ky for each control 
point and store the coefficients in the system analyzer. 

The deflection influence coefficients cy may be determined either by direct exper 
imental tests or by analytical methods. In this investigation an analytical approach was 
used. The analysis is based on linear shell theory of Sanders (ref. 11). A method of 
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solution using finite differences is outlined in reference 12 and a computer program using 
this approach was coded in reference 13. This program was chosen for this analysis 

f 

because it has general application and has been experimentally verified in reference 14. 
The program has the capability of accommodating 502 finite -difference stations along the 
shell meridian. Although the program is capable of including a large number of thermal 
control stations, only 10 stations were used, and each thermal load was permitted to span 
several node points in the finite -difference analysis. The midpoint of the thermal input 
may be selected as the control station for correction by the figure error sensor in the 
active control system. The temperature input necessary for the correction of any given 
deflection is, therefore, the amplitude of the same distribution used to determine the flex- 
ibility matrix. 

Using the computer program, the deformation at station i due to the application of 
a unit thermal input at station j was determined. The coefficients of this equation were 
determined as follows. Assume that the station of interest is number 1; then assign a unit 
amplitude (A in eqs. (1)) to Tj and zero to Tg, T3, . . ., TlO so that equation (2) 
reduces to 

5 1 = C 11 

6 2 = c 2 i 

63 = C 31 / ( 5 ) 


6 10 = e 101 J 


By applying a unit amplitude at each station separately and by calculating the resulting 
deflections with the computer program, each column of coefficients was determined. 


In order to adapt the computer program for use in this analysis, it was necessary 


to write subroutines describing the shell geometry and the thermal loading, 
eters necessary to describe the shell geometry are r, i c Oq, cog, and 

tions describing these parameters were formulated as follows: 


The par am - 

dcog 




Equa- 



+ f In 



( 6 ) 


(jOq = = =L= 

e 


(?) 
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( 8 ) 


_ 1 _ 4f 2 

5 R{ (4f2 +r 2 ) V2 

dfe, g _ _ 24rf 3 

d{ " («2 + r 2 

1ft* , 2t . (10) 

* r)/r 2 + 4f 2 

The program requires that the shell middle surface be defined by r = r(|). An explicit 
relation for r indicated in equation (10) cannot readily be written. Therefore, it was 
necessary to program r(£) by using a fourth-order Runge-Kutta approximation. The 
subroutine INPUT defining r(£) and the other shell parameters is given in appendix A. 

During the present investigation it was noted that some errors had been made in 
programing the thermal load terms presented in reference 13. These terms were in the 
subroutine FORCE. These terms were reprogramed and are included in the revised sub- 
routine FORCE given in appendix B. 



RESULTS AND DISCUSSION 


Deformation Sensitivity to Temperature Distribution 

Two patterns of control stations for the controlled temperature sources were inves- 
tigated. These patterns were for an equal meridional spacing and for an equal annular 
area spacing. The axial deflection at each of the 10 control stations for equal meridional 
spacing is shown in figure 5. These curves illustrate the deflection due to the application 
of a thermal input of unit amplitude at each station in turn as indicated in the figure. The 
deformations are a maximum at the inner unrestrained boundary (except for control sta- 
tion 7) and decrease uniformly to zero at the axially restrained outer boundary. It was 
noted previously that the tolerance criterion selected for this analysis was 0. 1 wavelength. 
At a wavelength of 0.5 pm, this mirror would have an allowable deformation 6 a of 
0.05 pm. The deformation created by a unit thermal input at station 10 is the largest 
and is nearly 600 times the allowable level. The axial deformation at each of the control 
stations for equal area spacing is shown in figure 6. These curves are similar to those 
shown for equal meridional spacing. With this spacing configuration, the maximum defor- 
mation created by a unit thermal input is larger than for meridional spacing and is over 
700 times the allowable level. 

The deflection influence coefficients were determined for both spacing configura- 
tions from the analytical data used to plot the curves shown in figures 5 and 6. These 
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coefficient matrices were inverted and the resulting stiffness matrices give the amplitude 
of the thermal input necessary to correct a given set of axial deflections at the control 
stations. In addition to the amplitude of the thermal input, the accuracy of the amplitude 
necessary for diffraction -limited operation may be determined. This accuracy is given 
by the product of the minimum contribution (minimum coefficient kjj) and the maximum 
allowable error, 6 a = 0.05 pm. An examination of the stiffness matrices indicates that 
the amplitude A of the thermal input must be controlled to within 1.0 mK for the case 
of equal meridional spacing and to within 0.23 mK for equal area spacing in order to 
insure diffraction limited operation. For either configuration the sensitivity of the defor - 
mations makes it necessary to employ a closed-loop control system to maintain the proper 
contour. Since the thermal inputs for the equal meridional spacing have lower control 
requirements, this spacing configuration is the preferred configuration of the two con- 
sidered herein. 


Constant Acceleration Loading 

To illustrate the use of the thermal active optics system the deformation due to a 
constant acceleration loading was examined. The axial deflection due to a 0.0 lg accel- 
eration load was determined by using the computer program of reference 13. The axial 
deformation as a function of radial location for this load is shown in figure 7. For this 
relatively light loading the maximum axial deflection exceeds the tolerance limit indi- 
cated in the figure by a factor of about 25. In order to correct this distortion it is nec- 
essary to introduce a deformation of equal magnitude and opposite sign by use of the 
thermal loads. By using the deformation of the O.Olg loading, the amplitude of the cor- 
rective thermal loads was determined from equation (4) for the case of equal meridional 



Figure 7 .- Axial deformation of ‘beryllium telescope mirror due to 
O.Olg constant acceleration load, f = 12.0 m. 



spacing. Then the computer program was used to determine the resulting displacement 
with both the acceleration load and the corrective thermal loads applied. The resulting 
deformation is shown in figure 8. Also indicated in the figure is the amplitude of the 
thermal input applied to each control point. The residual deformation is well within the 
tolerance limit; however, the deformation does not go to zero at all control points as 
might be expected. The residual deformation at the control points is due to errors intro- 
duced by the matrix inversion process and roundoff of the thermal inputs. It is significant 
that only low amplitude thermal inputs are necessary to correct deformations which exceed 
the tolerance limit by such a large amount. 


Amplitude of 

Control point thermal input, K 

1 .084 

2 .107 

3 . 092 

4 .095 

5 . 092 

6 .099 



Radial location, r/r 
' m 

Figure 8.- Axial deformation of beryllium telescope mirror subjected 
to O.Olg constant acceleration load and corrective thermal inputs, 
f = 12.0 m. 

CONCLUDING REMARKS 

An analysis for correcting axial deformations in a thin telescope mirror has been 
investigated. This analysis would be used in conjunction with a closed-loop control sys- 
tem to thermally correct surface deformations. The analysis has been demonstrated on 
a thin beryllium paraboloid with a focal length to diameter ratio of 4. Two patterns of 
control points for symmetrical thermal inputs were examined. To correct axisymmetric 
deformations the amplitude of the temperature distribution must be controlled to within 
1.0 mK for equal meridional spacing and to within 0.23 mK for equal area spacing. For 
either configuration the sensitivity of the deformations makes it necessary to employ a 
closed -loop control system to maintain the proper contour. Since the equal meridional 
spacing has a lower control requirement it is the preferred configuration of the two 
considered. 
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Deformations due to a constant acceleration load were examined and were found 
to exceed the allowable tolerance level by a factor of about 25. It was found that these 
deformations were readily corrected to well within the required accuracy by using rela- 
tively small amplitude thermal inputs. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 27, 1972. 
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APPENDIX A 


SHELL-GEOMETRY SUBROUTINES 


The subroutine INPUT necessary to calculate the shell geometry for the computer 
program of reference 13 is as follows: 

SUBROUTINE I NPUT (NMAX ) 

COMMON R ( 502 )/BL 1 /GAm ( 502 ) » OMT ( 502 ) . OMX I ( 502 ) * DEOMX ( 502 ) /BL3/NU . 

1 LAM, n. FALSI G. CHAR, DEL 

READ ( 5 « 1 )S»F,R(1 ) 

1 FORMAT (3E 10,5) 

DEL = f/ ( CHAR* FLO AT (NMAX-2 5 ) 

X=1 

A=CHAR 

NTLAS=NMAX— 1 

IE (R(l.) .EO. 0,0000 ) 2,? 

2 GAM ( i )=I »F32 1 
GO To 4 

3 GAM ( l ) = 1 • / ( R ( K ) *SQRT (R ( X ) **2*A**2/ (4 »*F**2 ) + l . ) ) 

4 OMT (i ) = 1 ./SORT (4. *F**2/A**2+R (K ) **2 ) 

OM X I ( 1 ) =4.*A*F**2/(4.*F**2+A**2*R (X 5**2 )**1 .5 
DEOMx ( 1 ) =-24 ,*A**3*R (X ) *F**3/ ( 4 . *F**2 + A**2*R < < > **2 ) **3 
CALL SURF AC ( 2 »f , 0 , 5*nFL 5 
DO 5 X = 3 , NT(_ AS 
CALL SUREAC (K ,F ,DFL 1 
=! CONTINUE 

CALL SURFAC (NMax.F,0,5*DEL ) 

DO 6 K = 2 » NMA X 

GAM (<)*! ./(R(X)*SQRT(R(K)**2*A**2/(4.*E**2)+1 • ') ) 

OMT (< ) = 1 ./SORT <4.*E**2/A**2+R(K )**2 ) 

OMXI (X ) =4.*A*F**2/ (4 ,*F**2 + A**2*R <K )**2 )**1 .5 
6 DEOMX (X ) =-?4.*A**3*R (X ) *F** 3/ ( 4 . *E**2+A**2*R ( X 5**2 5 **3 

rftupn 

END 
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APPENDIX A 


The subroutine SURFACE called by the subroutine INPUT is given as follows with 
the necessary function subprogram: 

SUBROUTINE SURFAC (K » F « DELTA ) 

COMMON R < 502 ) /BL3/NU ,L AM , N <EALS 1 G , CHAR , DEL 
REAL K1 tK2«K3»K4 
K1 =DELTA*VALUE (R (K-l ).F,CHAR) 

K2=DELTA*VALUE <R(K-1 )+0.5*Kl .F.CHAR) 

K3=DELTA*VALUE (R (K-I )+0.5*K2 ,F.CHAR) 

K4=DFLTA*VALUE (R (K-l ) +K3 1 F 1 CHAR ) 

R(K) =R (K-l ) + < 1 ./6. )* (K1 +2.*K2+2.*K3+K4 ) 

RETURN 

END 


FUNCTION VALUE (P» F « ChAR ) 
C=2.*F 

D=SQRT ( P*P*CHAR*CHAR+4 .*F*F ) 

VALUr=C/D 

RETURN 

END 
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APPENDIX B 


SUBROUTINE FORCE 


The corrected subroutine FORCE for the computer program of reference 13 is 
given as follows: 

subroutine: force (<» indb.nmax) 

c subroutine force this subroutine calculates the elements of the lower case 

C E-VECTOR AS DEFINED IN APPENDIX A OF REFERENCE ( 12 ) « AT THE STATION SPECIFIED 

c by the index k. 

common R <502 )/BLl/GAM (502 ) .0MT(502 ) «OMXI (502 ) »DEOMX(502 )/BL3/NU, 

1 L AM , N » EALS I G * CHAR ♦ DEL/BLA/CEE ( 4 ) 

REAL NU » LAM *N * L2 
REAL MSTP 
RA=R(K) 

GA=GAMCK) 

OX = OMXI (K > 

OT=OMT (K ) 

T = TFmP ( K » DEL ) 

BfMIS&ldftIKbL > 

DLT1 =DDELT (K » DEL ) 

PX 1 =PX ( K * DEL ) 

PT1=PT(K,DEL) 

P1=P(K,DEL) 

H=HHT(K.DEL) 

DH=DhHT(K.DEL) 

HRB=hRA (<»DEL ) 

DHRB = DHRA (K » DEL ) 

D 1 =1 ,-NU 

L?=LaM#*2 

EAL=pALSIG 

TSUBT=2 »*H*EAL/ (D1 *HRB )*T 

MSTP=CHAR*EALX ( 3 . *D1 )*( (t . 5/HRB-3 ./HRB**2+2. /HRB**3 ) * ( DLT 1 *H**3+3 * 
1*H**2*DH*DELT1 >+DELTi*H**3*DHRB/HRB**2* (-1 .5+6 ./HRB-6 • /HRB**2 ) ) 

CEE(4 )=CHAR*EAL*DELTl*H**3/(3.*Dl >*■.( 1 • 5/HRB-3. /HRB**2+2 • /HRB**3 ) 

CEE ( 1 ) --PXl +2.*EAL/' (Dl *HRB >* (H*DT+T*DH>-DHR3/HR3*TSUBT- 
1L2*D1 *GA*OX*CEE <4 ) 

CEE ( 2 ) = -PT 1 -N/RA*TSUPT-L2*D 1 *N/RA*OT*CEE (4 ) 

CEE ( 3 ) = -P 1 - ( OX+OT J*TSUBT-L2*D 1 *GA*MSTP+ 

1 L2*D 1 *CEE ( 4 ) * < OX*OT~ <N/RA)**2 ) 

RETURN 

END 
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